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We consider the possibility of adding noise to a quantum circuit to make it efficiently simulatable 
classically. In previous works this approach has been used to derive upper bounds to fault tolerance 
thresholds - usually by identifying a privileged resource, such as an entangling gate or a non-Clifford 
operation, and then deriving the noise levels required to make it 'unprivileged'. In this work we 
consider extensions of this approach where noise is added to Clifford gates too, and then 'commuted' 
around until it concentrates on attacking the non-Clifford resource. While commuting noise around 
is not always straightforward, we find that easy instances can be identified in popular fault tolerance 
proposals, thereby enabling sharper upper bounds to be derived in these cases. For instance we find 
that if we take Knill's [J] fault tolerance proposal together with the ability to prepare any possible 
state in the XY plane of the Bloch sphere, then no more than 3.69% error-per-gate noise is sufficient 
to make it classical, and 13.71 % of Knill's 7 noise model is sufficient. These bounds have been 
derived without noise being added to the decoding parts of the circuits. Introducing such noise in 
a toy example suggests that the present approach can be optimised further to yield tighter bounds. 



An important open problem in the field of quantum 
computation is to understand the effects of noise on com- 
putational power. In particular, for various noise models 
and various sets of universal quantum resources we would 
Ukc to determine the so called fault tolerance threshold Q , 
the level of noise that we can tolerate in our basic physical 
components before we lose the power to perform quan- 
tum computation. Most previous work on this question 
has focused on the problem of constructing lower bounds 
to such thresholds (see e.g. (H). This is usually achieved 
by constructing explicit fault tolerant procedures, and 
then estimating the level of noise that these schemes can 
tolerate. The construction of upper bounds on the other 
hand has received comparatively less attention, but has 
recently been the subject of an increasing number of in- 
vestigations. 

There are essentially two methods that have been used 
to derive upper bounds, which we may loosely describe 
as 'classical' and 'quantum' approaches. In 'classical' ap- 
proaches one tries to determine the noise levels required 
such that the noisy quantum computer can be efficiently 
simulated on a classical computer 

'quantum' approaches one tries to show that above a cer- 
tain level of noise the output is bounded away from ideal 
in some way @, 0, [13, U , l2| ■ A number of investigations 
along these lines have resulted in upper bound estimates 
relevant to various architectures and noise models un- 
der varying degrees of assumption and rigour. The two 
approaches could ultimately be of independent interest, 
as it is possible that for realistic noise models an 'in- 
termediate' form of computation could exist that is not 
as powerful as a quantum computer, but more powerful 
than classical. 

In this work we will follow the path of , and con- 



sider 'classical' threshold upper bounds for a specific class 
of fault tolerant quantum computational schemes - those 
built around a core set of Clifford resources 13| such 



as the CNOT, Pauli rotations, Hadamard, and prepa- 
ration/measurement in Pauli operator eigenbases. Any 
device consisting of such Clifford operations can be effi- 
ciently simulated classically, and so to perform quantum 
computation an additional resource is needed. Common 
choices include phase gates of the form 



Ui9) |O)(O|+expi0|l)(l|, 
or a supply of single qubit states such as 
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that are not eigenstates of Pauli operators. 

It is significant for our present discussion that devices 
built from Clifford operations alone can be efficiently sim- 
ulated classically. This is the content of the Gottesman- 
Knill theorem l3, The Gottesman-Knill theorem is 
interesting as in many other respects Clifford operations 
exhibit a great deal of non-classicality - for instance they 
can be used to demonstrate the GHZ paradox and may 
generate the long-range entanglement that is considered 
to be a precondition for efficient quantum computation 
[lit . Furthermore the entanglement that Clifford opera- 
tions generate is sufficient for quantum computation - to 
build a full quantum computer you only need to add the 
ability to perform single qubit operations to the Clifford 
set. 

The Gottesman-Knill theorem is important for this 
work as it can be used [1, [1| to construct classical thresh- 
old bounds using the following simple approach. Con- 
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sider a Clifford based architecture which is made univer- 
sal by the addition of a non-Clifford resource, which we 
refer to as R. If one could determine a noise level suffi- 
cient to turn R into an operation that lies in the convex 
hull of Clifford operations, then this noise level provides 
an upper bound to the fault tolerance threshold for those 
architectures, as such a noisy device can be efficiently 
simulated classically [22|. Although one would ideally 



like to derive threshold bounds that are valid for a wider 
variety of possible architectures than only those based 
upon Clifford operations, focussing on such architectures 
is nevertheless quite relevant, as most lower bounds have 
been derived for Clifford based schemes due to their sig- 
nificance for quantum error correction 14 1. 

Although a number of interesting threshold bounds can 
be derived with this approach, the noise models con- 
sidered in 0, are in fact be relatively weak because 
the Clifford operations are taken to be entirely noise 
free. This naturally leads to the question of whether 
the bounds can be improved by considering the (often 
more realistic) situation in which the Clifford operations 
are also subject to noise [l^. In this work we argue that 
this is indeed the case. By a straightforward modification 
- allowing the Clifford operations to be noisy, and then 
'commuting' the noise onto the non-Clifford operations - 
we argue that the previous bounds of 0, [1] can be im- 
proved much further for some important families of fault 
tolerant quantum computation. 

'Commuting' noise around a circuit is not always 
straightforward as when moved through entangling gates 
it can lead to unmanageable error correlations between 
different qubits. Hence in this article we analyse only 
comparatively easy instances where these correlations 
can be made to disappear or may be accounted for easily. 
Fortunately, however, it turns out that the tcleportation 
circuit is one such 'easy' instance, and so the approach 
works well for the various schemes (including the high- 
threshold proposal of Knill [l|) which use tcleportation 
as a primitive for 'injecting' non-Clifford states into the 
computation to make it universal. In such cases the prob- 
lem often reduces to simple geometrical considerations of 
the Bloch sphere. We expect that the results obtained 
here may be improved by developing more sophisticated 
techniques to tackle the correlations arising in more gen- 
eral situations. 

A summary of the upper bounds obtained in this pa- 
per is given in table (|T|. In some cases they are not 
too far from conjectured lower bounds for some proposed 
schemes. 



BASIC DEFINITIONS AND NOISE MODELS 

There are a variety of noise types and models that we 
will consider. In this section we define some of these 
noise models, and establish some notation. We will only 



consider computational models using qubits. 

• Clifford Operations. Consider the 'Clifford 
group' of unitaries, defined for any number of 
qubits, which is generated by CNOTs, Hadamards, 
and the 'S" gate S := |0)(0| Augment 
this group with measurements in the Pauli X, Y, Z 
eigenbases, and the ability to prepare any Pauli 
eigenstate. Any physical operation that can be 
generated by probabilistic application of these re- 
sources (where the probabilities can be efhciently 
sampled classically)will be called a Clifford opera- 
tion. 

• Phase gates. Phase gates, denoted U{6), are de- 
fined as: 



U{e) := |O)(O|+exp(z0)|l)(l| 



(3) 



In most of the paper we abuse notation and also 
use U to refer to the quantum operation (on density 
matrices) that U corresponds to. 

• Phase states. Phase states, denoted \0) are de- 
fined as: 
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— (|O)+exp(z0)|l)). 



(4) 



They lie in the plane of the Bloch sphere containing 
the Pauli X, Y eigenstates. 

• Probabilistic Application of a transforma- 
tion Q. Let Q be a quantum operation acting on 
a quantum system. Q could be a unitary or a com- 
pletely positive (CP) map. We define to be the 
quantum operation "with probability (1 — t) apply 
the identity operation to the system and with prob- 
ability t apply the Q operation" , i.e. 



(5) 



Usually we will consider Q to be a Pauli operation, 
or some form of depolarising operation. 

• Opposite Noise. Consider a known, given, qubit 
state a. Define to be the pure state in the polar 
opposite direction to a in the Bloch sphere. We 
define Of to be the single qubit operation "with 
probability {1 — t) apply the identity operation to 
the system and with probability t throw it away 
and replace it with the pure state in the opposite 
direction to cr in the Bloch sphere" , i.e. 



(6) 



This operation will be used to construct adversarial 
noise upper bounds. Note that it is a valid quantum 
operation as a is known and fixed, and does not 
depend upon p. 
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• Knill's noise model. The noise model defined in 
the paper by Knill [ij is parameterized by a noise 
strength 7 in the following way: a preparation of 
a qubit state ip prepares the orthogonal pure state 
with probability 47/I5, a measurement of the Pauli 
Z operator is preceded by a Pauli X rotation with 
probability 47/I5, a measurement of the Pauli X 
operator is preceded by a Pauli Z rotation with 
probability 47/ 15, a single qubit unitary is followed 
by one of the Pauli X, Y, Z rotations each with 
probability 47/I5, a CNOT is followed by one of the 
15 non-Identity Pauh products {1(g) X, Y®I, Z®X 
etc) each with probability 7/ 15. 

• Independent Depolarising noise, with 
strength t. After every non-trivial unitary oper- 
ation, and before each non-trivial measurement, 
apply the noisy operation iV/^, where D is the 
totally depolarising operation: 

D{p)^\{p + XpX + YpY + ZpZ) (7) 

• Simultaneous Depolarising noise, with 

strength t. After each non-trivial gate acting 
on k qubits, each possible non-Identity Pauli 
product is applied to those qubits with probability 
<:/(4" — 1). In particular this means that after 
a CNOT, for example, both qubits are 'jointly 
depolarised' rather than independently on each 



EPG. Probabilistic Error-per-gate Noise, 
with strength t. This is the noise model con- 
sidered in many of the more stringent and rigor- 



Each 



ous lower bound estimates, such as 
operation (including memory steps required while 
waiting for measurement outcomes elsewhere) fails 
with probability i, but the precise manner of fail- 
ure is not specified other than being probabilistic. 
It can in principle be adversarial. In particular the 
output of a /c-qubit gate can undergo some cor- 
related, joint, error. However, correlations in the 
noise affecting non-interacting qubits arc not al- 
lowed in this noise model. 



THE BASIC IDEA 

To illustrate the approach it is useful to begin with 
a specific example analyzed in [^. In that paper it was 
shown that if the Clifford operations are augmented with 
the ability to do single qubit phase gates and if each 
non-Clifford gate U{9) in the circuit is immediately fol- 
lowed by an unwanted Z rotation with probability p, then 
roughly p ^ 14.6% of noise is sufficient to make the cir- 
cuit classically tractable (this is regardless of which phase 



gates are used). However, let us now consider how this re- 
sult may be applied to more realistic noise models where 
Clifford operations are also subject to noise. As an initial 
case let us assume the following noise model: each qubit 
undergoes a dephasing operation of the form (see eq. 
([5])) after every non-trivial gate, after every qubit state 
initialisation, and before every measurement. 

Consider one particular qubit wire in a particular fault 
tolerant circuit. If time increases from left to right 
the qubit might be about to enter a sequence of (ide- 
alised) operations as follows: 



CUUC 



(8) 



where C denotes a Clifford operation (which may involve 
more than one qubit, and could be a state preparation), 
and U denotes a non-Clifford phase gate. The real noisy 
version of this wire will be: 



C TVf U TVf U TVf C TVf 



(9) 



By a little circuit manipulation we can reexpress this cir- 
cuit as one in which the noise is concentrated on the 
non-Clifford gates. The first step is to be generous and 
remove the noise acting in between the two Us, and on 
the last Clifford gate, to give (as is itself in the con- 
vex hull of Clifford operations, dropping any instances of 
it will not affect the validity of our arguments): 



CN^UU TVf C. 



(10) 



Now we can replace UU with [/', as the composition of 
two phase rotations will just be another phase rotation. 
This gives 



cn^u'n; 



C. 



(11) 



Hence we see that by dropping some of the noise we can 
reduce a circuit of Clifford operations and phase rotations 
to a circuit where (apart from the noise) every second 
gate on any individual wire is a Clifford operation. Now 
as our noise operation iV^ commutes with the t/, this 
equation reduces to: 



CU' {X^fC. 



(12) 



This is a circuit where the Clifford operations are per- 
fect, but we now have two lots of noise attacking the V 
non-Clifford gate. In some cases we will only have one 
U between two Clifford operations, and in some cases we 
may have many [/s of different types. We may also have 
a I] immediately after the preparation of a fresh qubit, 
or immediately before a measurement. However, in all 
these cases the reduction still works. In general it is not 
possible to argue that there are more than two lots of 
noise attacking each non-Clifford operation, as in princi- 
ple every second operation on every qubit wire could be 
a non-Clifford operation. 
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Hence we may now ask how high p must be for two lots 
of noise to take an arbitrary phase rotation U' into the 
set of Chfford operations - this is in contrast to the cal- 
culations of [3| which used only one lot of noise. Solving 
this problem is a trivial extension of the methods used in 
\4i. One must solve 



m 



(1 - 2pf = (1 - 2q) 



(13) 



for p, where q is the minimal noise level required for a 
single application of the noise. Using the precise value 
for q found in [5| we find 



1 f 1 

Solving for p hence gives 



P 




14.64%. 



7.955%. 



(14) 



(15) 



So this result may be summarized as follows: for an ar- 
chitecture based upon Clifford operations, phase gates, 
and noise where each qubit undergoes Np independently 
after every non-trivial gate acting upon it (or before a 
measurement), then a noise level of 7.955% is sufficient to 
make the circuit tractable classically. Some authors pre- 
fer to use a different parameterization of the dephasing 
noise model, where instead of applying Np , one applies 
the totally dephasing operation p — > l/2{p + ZpZ) with 
probability p. The two cases are trivially related, and so 
the above bound can be reexpressed for this noise model 
as 




15.9%. 



(16) 



These improvements on the bounds of Q have been possi- 
ble because noise from a neighbouring Clifford operation 
can be shifted onto the non-Clifford resource. It is easy 
to see that this attack could lead to improvements for 
any Clifford based architecture provided that the noise 
can be shifted around the circuit easily. In the example 
considered above, this 'shifting' was possible because the 
Np noise commutes with any U{9). Of course in general 
the noise models won't commute with the non-Clifford 
resource so straightforwardly, but as we shall see later, 
there are important examples of fault tolerance schemes 
in the literature for which this approach can be applied 
quite effectively. 
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FIG. 1: A possible scenario in a fault tolerant circuit. 

different physical qubit wires. However, a similar ap- 
proach can be used to derive a bound valid for an adver- 
sarial error-per-gate model, where each non-trivial op- 
eration is subject to probabilistic adversarial noise, but 
different physical qubits undergoing the same multiqubit 
gate (e.g. a CNOT) can have correlations in noise at 
that location in the circuit. The adversarial error-per- 
gate model is more commonly adopted in estimates of 
rigorous lower bounds, and hence is an important model 
to consider. 

In this case the analysis must be changed slightly. A 
circuit such as that in figure ([T]) may occur, where the 
same Clifford operation, in this case the CNOT, touches 
two non-Clifford operations. In this case at locations a, b 
wc choose to apply the noisy operation N^ ® Nf where 
t is chosen so that the overall probability of any error at 
locations a,b (of figure ^) together, + 2t(l — t), is 
equal to p, i.e. 



p = f + 2i(i ~t)=^t = i- 



P- 



(17) 



Now the noise affecting each non-Clifford operation in 
circuit [T] is effectively independent. As in this section we 
assume the non-Clifford operations to be a phase rotation 
gate (a similar argument follows through for preparation 
of a phase state, with precisely the same numerical val- 
ues), our goal is to work out when 

1 



(l-2p)(l-2(l- v/r^)) = l-2g=-^ (18) 



which has the solution for p of: 

p = 10.41008383%. 



(19) 



So under an adversarial EPG model, any circuit con- 
sisting of noisy Clifford operations and phase states and 
phase gates can have a fault tolerance threshold no higher 
than 10.41008383%. 



EPG BOUNDS FOR CLIFFORD OPERATION 
AND PHASE GATE OR PHASE STATE 
SCHEMES 



INDEPENDENT DEPOLARISING THRESHOLDS 
FOR GENERAL NON-CLIFFORD UNITARIES 



The bounds of the previous section have been derived 
for an adversarial noise model which is independent on 



Depolarising noise commutes with all single qubit 
unitarics, so provided that the non-Clifford operations 



are single qubit unitaries the above reasoning applies 
straightforwardly. Although in this section we now al- 
low U to be any non-Clifford single qubit unitary, let us 
in particular take it to the most robust gate to depolaris- 
ing noise, which in is shown to be the tt/S phase gate 
anyway. Let us ask what p is required such that (where 
as before time increases from left to right, / denotes the 
identity operation, and U represents the quantum oper- 
ation acting by conjugation on density matrices, not just 
the unitary matrix itself) 



t/(V4)«)^ 



U{TT/4){{l-p)I + pDf 
{l-pfUi7r/4) + il-{l 



pf)D 



is a Clifford operation (the second line of this equation 
follows from the fact that D'^ = DU = UD = D as D 
is the depolarising operation). We may now apply the 
bound obtained in J5|] (a related result has been derived 
independently by [8|) to give that the minimal such p 
satisfies: 



{i-pr = i- 



6 - 2V2 



= 100% - 45.3% 



hence 



p=l 



6-2V2 
7 



26.05% 



(20) 



(21) 



This upper bound applies regardless of the single qubit 
non-Clifford gates available and the fault tolerance meth- 
ods used, provided that the basic gates are Clifford + 
single qubit unitaries. 

Ben Rcichardt [lij has independently previously per- 
formed a related analysis for simultaneous depolarising 
noise. The calculations arc more difficult because in that 
case the noise docs not commute straightforwardly. How- 
ever, it is expected that an analysis of the simultaneous 
case would lead to upper bounds of a similar magnitude 



SPECIALISING TO TELEPORTATION 
STATE-INJECTION SCHEMES 

If general bounds are required for architectures involv- 
ing Clifford and non-Clifford operations, then it is gener- 
ally not possible to shift more than two lots of noise onto 
each non-Clifford operation. This is because in principle 
every second operation in a given qubit wire could be a 
non-Clifford operation. 

However, in many important proposals for fault tol- 
erant quantum computation the non-Clifford operations 
are always surrounded by specific configurations of Clif- 
ford operations. This is because non-Clifford operations 
are often introduced into the computation using very spe- 
cific ancilla preparation constructions. For example, one 
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FIG. 2: These circuits show two possible state injection set- 
tings. The boxes denote gates, the triangles are either prepa- 
ration (of the -1-1 eigenstate) or measurement in the X,Z bases, 
and the collection of wires at the bottom left of each circuit 
represents the postselected decoding process. The letters 'nc' 
denote 'non- Clifford'. The top circuit is for scenarios involv- 
ing single qubit preparation as the non-Clifford resource, the 
second is for single qubit gates. Each diagram has locations 
numbered from 1 to 6, with the lower figure having an extra 
location 7. In the text we use these numbers to define noise 
models for both the upper and lower circuit in this figure. For 
state injection protocols these circuits represent the only way 
in which non-Clifford operations enter the computation. 



method that we will consider here is the technique of 
'state-injection' - see figure ([2|). This method involves 
firstly creating a physical non-Clifford qubit, either by 
direct access to a source of such qubits, or by applying a 
gate such as the tt/S gate to a suitable Pauli-eigenstate. 
This qubit is then 'teleported' into an error correcting 
code, by first decoding one half of an encoded bell pair to 
the physical level (post-selecting on no errors), and then 
performing ordinary teleportation using CNOTs and X,Z 
measurements. Because the teleportation circuit imme- 
diately around the non-Clifford resource has around 5 to 
7 non-trivial error locations (depending upon the precise 
model), one can shift upto 7 lots of noise onto the non- 
Clifford resource, and obtain much lower upper bounds. 
By allowing noise at these locations, and then remod- 
elling this as effective noise acting only upon the non- 
Clifford resource, it is possible to strengthen the bounds 
derived in the previous sections. 

To illustrate the approach, let use consider what hap- 
pens when in the top circuit of fig.® we apply various 
types of Pauli error at locations 2 to 6. The following ar- 
guments give some rules for 'shifting' these Pauli errors 
to location 1: 

1. Application of Z at location 6, and no errors 
elsev^rhere. Because location 6 is immediately fol- 
lowed by a Z measurement, this case is essentially 
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equivalent to no error as it is 'absorbed' by the 
measurement. 

2. Application of X at location 6, and no er- 
rors elsewhere. An X on the target wire com- 
mutes with the CNOT, and so the X can in fact be 
commuted through to location 4. 

3. Application of Pauli error at location 4, and 
no errors elsewhere. Now because the CNOT 
and subsequent X, Z measurements implement a 
Bell measurement, and because each projector B 
onto a Bell state satisfies (P ® /)i?(P ® /) = (/ ® 
P)B{I ® P) for any Pauli operation P, a Pauli P 
at location 4 is equivalent to a P acting at location 
1. Putting this together with the previous point 
shows that an X at location 6 is equivalent to an 
X at location 1. 

4. Application of Y at location 6, and no errors 
elsewhere. A Y error as a quantum operation (not 
as a matrix) is equivalent to a Z error and an X 
error. The Z can be absorbed by the Z measure- 
ment on the lower wire, leaving an X error which 
can be moved to location 1 by the previous points. 
Hence a F at location 6 is equivalent to an X at 
location 1. 

5. Any Pauli errors at locations 4,5,6. Because 
Pauli matrices either commute or anticommute, 
when viewed as quantum operations Pauli errors 
actually commute with each other - consider for 
example the identity XZpZX = ZXpXZ . Hence 
the previous rules can be applied to any set of 
Pauli errors acting at locations 5,6 - effectively re- 
expressing them as Pauli errors acting with various 
probabilities at location 1. 

6. Any Pauli errors at locations 2,3. Any X er- 
rors can be absorbed by the X measurement on the 
top wire, any Z errors can be commuted through 
the CNOT to location 1. 

We consider two variants of the state injection schemes. 
In state resource variants, the non-Clifford resource is a 
pure qubit state as in the top circuit of figure (l2|), in 
gate resource variant, the non-Clifford resource is a single 
qubit unitary, as in the bottom circuit of figure (|2|). 

The above rules can be applied to any configuration 
of Pauli noise at locations 2-6 to shift it all to location 
1, where it can attack the non-Clifford resource. As the 
non-Clifford resource at location 1 is effectively a state 
in the Bloch sphere, we can solve the relatively simple 
problem of how much of this noise forces the state to 
enter the Clifford 'octahedron' (c.f. [^, [l3|) formed from 
the convex hull of Pauli eigenstates. In the case of upper 
bounds for an EPG model, we are also free to try to pick 



noise 




FIG. 3: A 2D projection of the Bloch sphere. The goal of the 
calculations in this paper is to determine the noise strength 
required for an input state to enter the octahedron formed 
from the convex hull of the Pauli eigenstates. 

the most adversarial noise we can. Consider a pure Bloch 
vector 

V j (22) 

in the positive octant of the Bloch sphere (i.e. x,y,z > 
0). A Z error flips the sign of x,y, an X error flips the 
sign of y, z. Our goal is the find the minimal noise, for a 
given noise model, which takes the input Bloch vector to 
an output one 









Noise : 




y' 









satisfying x' +y' + z' = 1, which is the equation of the face 
of the octahedron in the positive octant (see e.g. figure 
©). 

BOUNDS FOR THE KNILL NOISE MODEL AND 

TELEPORTATION STATE-INJECTION, FOR 
PHASE STATE AND PHASE GATE RESOURCES. 

To be perfectly clear, it is important to specify pre- 
cisely how we apply Knill's noise model to the circuits in 
figure [21 In the top circuit of figure [2] we apply: a Z at 
locations 1,3 with probability 47/15; an X at location 6 
with probability 47/15; and at locations 2,5 considered 
together we apply each of the 15 non-identity pairs of op- 
erators chosen from /, X, Y, Z , each pair with probability 
7/15. Note that we have kept location 4 error free. The 
noise at location 4 will be determined by the decoding 
circuit that feeds it, and so in order to obtain general 
bounds independent of the codes used we are adopting a 
noise model that is strictly weaker than Knill's. Later we 
will discuss the effect that noise in the decoding circuits 
might have. 
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In the bottom circuit of figure [2l on the other hand, 
we apply the noise model: at 3,7 we apply a Z with 
probability 47/15; at 6 an X with probability 47/15; at 1 
we apply a X,Y and Z each with probability 47/15; and 
at locations 2,5 considered together each of the 15 non- 
identity pairs of operators chosen from /, X, Y, Z, each 
pair with probability 7/15. Again we keep location 4 
error free. 

Using this noise model leads to the following effective 
transformation of the input Bloch vector: 



X' 

y' 



= l- 



I67 
T5" 



1 - 



87 
T5 



(1-87/15)0; 
(l-87/15)2y 
(1 - 87/15)2 



where n ~ 1 for the top circuit of figure ([2]), and n = 
2 for the lower circuit of figure Our goal is hence 
to determine, for a given input resource, the minimal 7 
such that the output Bloch vector lies on the face of the 
octahedron. 

For example, if we assume that the ideal state entering 
location 1 in both circuits is |7r/4), either because that is 
the state prepared, or because the non-Clifford unitary 
is the 7r/8 gate, then find the solution: 

(1 - 167/15)(1 - 87/15)2(1 + (1 - 87/15)) ^ V2 

7 - 13.6861%. (24) 

for the top circuit of figure [21 and 

(1 - 167/15)2(1 - 87/15)2(1 + 87/15)) = V2 

7 - 9.5858%. (25) 

for the lower circuit of figure [2l Similar equations can 
be derived and solved for any possible phase gate and 
phase state resource. It is not difficult to solve for the 
minimal 7 that leads to an output vector on the face 
of the octahedron. Numerically scanning through these 
solutions suggests that the |7r/4) and C/(7r/4) resources 
are not actually the most robust non-Clifford phase state 
or phase gate resources in this setting, although they are 
very close. If we allow all phase state and phase gate 
resources respectively, the upper bounds become: 
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for phase states, and 



7 ' 



13.71% 



9.59% 



(26) 



(27) 



for phase gates. Later we present the values obtained if 
any single qubit state or gate is permitted as the non- 
Clifford resource. 



BOUNDS FOR AN EPG NOISE MODEL AND 
FOR PHASE GATE OR PHASE STATE 
TELEPORTATION STATE INJECTION 



ble, yet is sufficiently simple to analyse. It is quite likely 
that the noise that we pick at each location is not the 
most adversarial (in the sense of pulling the teleporta- 
tion circuit into a Clifford operation) within the EPG 
constraint, but this would require further analysis. 

We choose the following noise for the top circuit of 
figure ([2]). Assume that the input non-Clifford state is 



Apply Oj, 



at location 1; at locations 3,4; 

mm 



X 



at location 6; and replace the CNOT with Op (see 
eq. (O) on the top wire followed by a CNOT. This gives 
the following transformation on the input Bloch vector 
(which has z = as it is a phase state): 



(;)-(l-2p)^((l-p)(l-2p)-p)((^_^2p). 



For an input state \9 ~ 7r/4) , the minimal p such that the 
output Bloch vector lies on the face of the octahedron is: 



p = 3.68124%. 



(28) 



Numerics again suggest that \6 ~ 7i'/4) is not the most 
robust state in this setting, and in fact 



p = 3.69% 



(29) 



is a noise level sufficient to take all possible phase states 
into the octahedron. 

For the lower circuit of figure ^ we choose the fol- 
lowing noise. Assume that the input state is |-|-), the 
-f-l eigenstate of A", and that U is the non-Clifford phase 
gate. Apply oj,^^^'*'' at location 7; O^'"''^^^''^ at location 
1; at locations 3,4; iV^ at location 6; and replace 

the CNOT with Op on the top wire followed by 
a CNOT. This gives the following transformation on the 
output Bloch vector (which has z = as it is a phase 
state): 



(l-2p)2((l-p)(l-4p + 2p2)_p)f _^ 



(1 - 2p)y 



For an input state U{9 ~ 7i'/4) the minimal p such that 
the output Bloch vector lies on the face of the octahedron 
is: 



p = 3.00339%. 



(30) 



Numerics suggest that C/(6' = 7r/4) is not the most robust 
phase gate in this setting, and in fact 



p = 3.01% 



(31) 



To obtain upper bounds for an adversarial EPG model, 
we will pick noise which appears as detrimental as possi- 



is a noise level sufficient to take all possible phase states 
into the octahedron. 
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BOUNDS FOR TELEPORTATION INJECTION 
WITH GENERAL NON-CLIFFORD STATES AND 
GATES 

It is not difficult to perform the analysis of the previous 
two sections using general single qubit non-Clifford gate 
and state resources, and then to numerically calculate 
the most robust of these resources. For the Knill noise 
model we find that for general unitary gates 

7 = 15.19% (32) 

is sufficient to turn the teleportation state-injection into 
a Clifford circuit, whereas for general states 

7 = 21.78% (33) 

is sufficient. The most robust non-Clifford state entering 
the circuit appears to be close to, but not exactly the 
same as, the so-called |T) state. 

For an EPG model we use the same noise as for the 
phase gates/states in the previous section, except at loca- 
tion 4 instead of applying we apply A''^. Numerical 
analysis of the resulting equations finds that: 

p = 5.03% (34) 

is sufficient to turn the teleportation state-injection into 
a Clifford circuit, whereas for general states 

p = 6.31% (35) 

is sufficient. 



POTENTIAL EFFECTS OF THE DECODING 
CIRCUITS. 

So far, wc have considered solely the injection part of 
the circuit as depicted in fig. [2l and have not attempted 
to understand the effects of noise in the decoding circuit. 
Indeed, within the framework of Knill's noise model wc 
have not even allowed noise at location 4. The precise 
form of the noise generated by the decoding circuit is 
difficult to determine due to the complex structure of 
the encoding and decoding networks, which may include 
many steps of concatenation. This brings with it two 
problems. Firstly, the entangling gates in the decoder 
typically generate correlations in noise which mean that 
the problem becomes a multiqubit problem, rather than 
a simple geometrical problem on one or two qubits. Sec- 
ondly, the postselection steps have an effect on the noise 
profile, and this is difficult to calculate. 

However, we believe that a careful analysis of the de- 
coding circuits will lead to improvements in the bounds 
that we have presented. In this section we present some 
very rough indications of the level of improvement that 




(2)) = Depolarizing error (G) = CNOT error 

@ = Measurement error, either Z or X depending on subsequent measurement 
= Preparation error as Z error 



FIG. 4: This circuit shows the state-injection circuit of fig. 
[2] (lines 1 and 2) but now we explicitly include the action of 
part of a decoding circuit with the Knill error model. The 
input of the circuit is the state (|i/)o)|0) + |'!/'i)|l))/\/2 where 
the state will be decoded by an error free circuit to \i) in 
the second qubit. Errors from previous stages of encoding and 
manipulation are modelled by depolarizing errors. The error 
free state-injection would then lead to the preparation of the 
state \'4>) in the Output qubit at the bottom. Determining 
the minimal error rate that ensures that the Output is a sta- 
bilizer state provides an upper bound on the fault tolerance 
threshold. 



might be expected, although a more careful analysis is 
left to another occasion. 

To get a feeling for the effect of the errors that are in- 
troduced by the decoding circuit that immediately pre- 
cedes position 4 we have carried out an analysis of the 
circuit in fig. |31 which corresponds to one level of de- 
coding prior to the injection circuit for a particular code. 
The input of the top arm of the teleporter is a phase 
state proportional to |0) + cxp(m/4)|l). This circuit is 
a simplified version of Fig. 9 of [20|. The simplification 
consists of neglecting essentially all errors that may occur 
in the preparation and encoding, as well as the decoding 
of the lower half of that circuit of Fig. 9 of [2^. We 
also model the error on the incoming qubits in fig. [¥]via 
depolarizing errors (it is in this assumption that most of 
the complexity is buried). 

The computation of the error threshold in this ap- 
proach is possible analytically because of the small cir- 
cuit size. With the help of a computational mathematical 
package we find the threshold to be the sole real root of 
the polynomial 

2 /2(ec) 
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TABLE I: A summary of upper bounds presented in this work. 
'Injection' refers to teleportation state-injection. 'NC means 
non-Clifford. 



NC resources 


Method 


Noise Model 


Upper bound 


phase gates/states 


any 


independent Np 


7.96 % 


phase gates/states 


any 


EPG 


10.41% 


all gates 


any 


indep. depolarizing 


26.05% 


phase gates 


injection 


KniU 


9.59% 


phase states 


injection 


KniU 


13.71% 


phase gates 


injection 


EPG 


3.01% 


phase states 


injection 


EPG 


3.69% 


all gates 


injection 


KniU 


15.19% 


all states 


injection 


KniU 


21.78% 


all gates 


injection 


EPG 


5.03% 


all states 


injection 


EPG 


6.31% 



with 

/i(ec) 

/2(ec) 



= (16ec- 15)^(225- 180ec + 326^)2 X 

X (16875 - 40500ec + 468006^ - 23040ei? -f 4096e^) 

= 1125\/2(576650390625- 3536789062500ec 



+ 11768793750000e^ 
+ 32367600000000e^ 
+ 18141511680000e^ 
+ 1887436800000e^ - 
+ 17179869184e;^") 



- 24002325000000e;! 

- 29499033600000e^ 

- 7375159296000e^ 
273804165120e? 



and wc find 



threshold 



9.2888% 



(36) 



One can try to understand how critical is the choice of 
depolarising errors entering the decoding steps. To es- 
timate this we also analysed what happens if the depo- 
larising errors in the central four wires of figure ^ are 
removed. In this case the bound becomes 



threshold 



10.0638% 



(37) 



Both these bounds are quite close to the bound that 
would be obtained by instead modelling the noise from 
the decoding circuit as a depolarising error at location 4 
(for which the bound would be 9.59%). All these num- 
bers should be compared to 13.69%, which is the number 
obtained for precisely the same phase state input with no 
noise whatsoever outside locations 1-3 and 5-7 of Fig. [51 
A more rigorous analysis of the decoding circuits re- 
quires much more effort and may be reported elsewhere. 
However, these rough calculations give some indication of 
the improvements that might be expected in the upper 
bounds that we have obtained so far. 



DISCUSSION, CAVEATS, AND CONCLUSION 

We have discussed 'attacks' on fault tolerant quantum 
circuits involving Clifford operations and extra resources, 
with the intention of adding the smallest amount of noise 
possible to make the circuits efficiently simulatable classi- 
cally. The approach is simple - to shift noise from neigh- 
bouring Clifford gates onto the non-Clifford resources. 
The approach works best for situations where this 'shift- 
ing' process can be done easily, as happens for teleporta- 
tion based state injection schemes. There are many other 
fault tolerance proposals in the literature (e.g. meth- 
ods built around cluster states [2l[) that involve a few 
non-Clifford resources surrounded by many Clifford op- 
erations - and so in such cases our approach could also 
provide improved upper bounds, depending upon the pre- 
cise manner in which the 'non-Cliffordness' is 'injected' 
into the rest of the circuit. 

We are certain that the upper bounds obtained in this 
work can be optimised further, particularly through con- 
sideration of the decoding circuits. The rigorous part 
of the analysis performed in this work does not consider 
the decoding circuit present in state-injection schemes, 
except for a noise step at location 4 in the error-per-gate 
noise model. The analysis of the teleportation part of the 
state-injection is straightforward because the teleporta- 
tion circuit simplifies correlated Pauli errors, reducing 
the problem to a single-qubit one. A more sophisticated 
attack that also apphes noise to the decoding steps should 
lead to tighter bounds, especially for Knill's noise model. 

It is also important to note that teleportation is not 
the only method of state injection. If the non-Clifford 
resource is either a measurement or a unitary, a non- 
Clifford measurement may be implemented directly on 
the decoded wire. This method involves far fewer error 
locations outside the decoder, and hence is less suscepti- 
ble to our method of attack. Hence the bounds derived 
in the lower part of table (|T| are not generally applica- 
ble to all stabilizer schemes, but are intended more as 
an approach that can be fairly effective for specific fault 
tolerance proposals. 

A summary of the upper bounds that we have obtained 
is given in table (jT|. Although the bounds that we have 
derived are lower than previous rigorously established up- 
per bounds, our bounds are not typically as general as 
most previous results make fewer assumptions about the 
architecture (and often involve incomparable noise mod- 
els). For comparable fault tolerance schemes the esti- 
mates conjectured in [7] are lower than the values that 
we have obtained, and future work will be needed to see 
whether further optimisation of our approach will go as 
low. However, the approach presented here is relatively 
simple, makes no assumptions (other than assuming that 
quantum computation is not classically tractable), and 
has a slightly different aim as it deals with 'classical' 
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bounds rather than 'quantum' ones. 
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